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Abstract 

The p spectral function at finite temperature calculated using the real-time formalism of thermal field 
theory is used to evaluate the low mass dilepton spectra. The analytic structure of the p propagator is 
studied and contributions to the dilepton yield in the region below the bare p peak from the different cuts 
in the spectral function are discussed. The space-time integrated yield shows significant enhancement in the 
region below the bare p peak in the invariant mass spectra. It is argued that the variation of the inverse 
slope of the transverse mass (My) distribution can be used as an efficient tool to predict the presence of two 
different phases of the matter during the evolution of the system. Sensitivity of the effective temperature 
' obtained from the slopes of the Mj spectra to the medium effects are studied. 

> 
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^ 1 Introduction 

^ \ The latest results from relativistic Au + Au collisions at Relativistic Heavy Ion Collider (RHIC) have indicated 
• that the matter created in the initial stages in relativistic heavy ion collisions (HIC) might be in the form 
. of quark gluon plasma pQ. This is supported by the observation of high transverse momentum (pr) hadron 
T 7"! \ suppression (jet-quenching) in central collisions compared to the binary-scaled hadron-hadron interactions [2]. 
# ^ ' The observation of large elliptic flow have also indicated the possibility of rapid thermalisation of such high 
^ . density matter [3]. 

Although the small production cross-section of lepton pairs leads to a lower yield, the electromagnetic probes 
are well suited to probe the local properties of the transient form of matter produced in nuclear collisions at 
ultra-relativistic energies as they leave the system almost unscathed. They emanate at all stages and thus are 
expected to map the temperature profile of the evolution. Because of the backgrounds coming from different 
stages and the low yield, extraction of signal from the background becomes a daunting task [4]. Theoretically, 
the prediction of the yield depends on the evaluation of the production rate as well as the scenario of space-time 
evolution that one employs. The rates of production from the QGP is controlled by QCD but those from 
hadronic matter depends on the hadronic interactions that one considers. In the low mass region, the rate of 
dilepton production is controlled by the spectral functions of the vector mesons, specially the p and hence the 
modification of the p spectral function determines the yield of lepton pairs in this region of invariant mass. 

A number of authors have analysed the dilepton spectra from heavy ion collisions; the treatments differing 
both in the construction of the p spectral function as well as the space time evolution scenario employed. This 
includes the nature of phase transition, the equation of state as well as numerical values of the parameters like 
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the initial temperature, the thermalisation time, the phase transition temperature as well as the chemical and 
kinetic freeze-out temperatures. We do not attempt to review or summarise the considerable amount of work 
which has been done on this topic except to mention the most recent few. The NA60 experiment at the CERN 
SPS measured dimuon pairs in In-In collisions in which an excess was observed over the contribution from 
hadronic decays at freeze-out in the mass region below the p peak [JJ. This was attributed to the broadening 
of the p in hot and dense medium [5], in contrast to the earlier data from the CERES collaboration [6] which 
is unable to differentiate between the broadening and the pole shift of the p spectral function [7J. The NA60 
data for the entire (measured) invariant mass range is reproduced by taking into account dilepton productions 
from Drell-Yan processes, qq annihilation, thermally broadened in-medium p, decays of p at the freeze-out 
surface and primordial p produced from the initial hard scattering [$]. The dependence of thermal emission 
on transition and chemical freeze-out temperatures is also highlighted there. The dilepton yield evaluated with 
the in-medium spectral functions of p and u mesons deduced from empirical forward scattering amplitudes [9] 
does not reproduce the data well at the low invariant mass (M < 0.5 GeV) region. The PHENIX experiment 
reported a substantial excess of electron pairs in the same region of invariant mass [10J. The data has been 
investigated by several groups e.g. 112j . The yield in all these cases have remained insufficient to explain 
the PHENIX data |13j . Thus the issue of low mass lepton pair yield still remains an unsettled issue. 

In the literature the modification of the p propagator comes from (a) the thermal modification of the decay 
width into pion pairs and (b) collisional broadening due to scattering with the excitations in the thermal 
medium [5]. In this work we have studied dilepton production at RHIC and LHC (Large Hadron Collider) 
energies using a spectral function evaluated in the real time formulation of thermal field theory with the 
interaction vertices taken from chiral perturbation theory which is the low energy effective theory of QCD. 
From the discontinuities of the p self energy associated with the branch cuts in the complex energy plane we 
provide a unified description of the apparently different scattering and decay processes in the medium [14] (for 
a different approach see [15\ [5] ) . The resulting spectral function of the p at non-zero three- momentum shows 
significant broadening with no appreciable change in the pole position. Since the dilepton spectra is proportional 
to the vector meson spectral function, the p in particular, we will attempt to bring out the medium modifications 
through different aspects of the dilepton spectra (see [5j [16] for review). 

For the space-time evolution we have used ideal relativistic hydrodynamics. There are quite a few parameters 
which go as inputs into this scheme e.g. the thermalisation time or the initial temperature. Though studies on 
elliptic flow of matter produced in such collisions indicate a rapid thermalisation [T7] , analyses using second order 
transport coefficients with conformal symmetry indicate sizable uncertainties [18] in its determination. Then 
comes the issue of transition temperature in addition to the chemical and kinetic freeze-out temperatures. The 
equation of state, in particular the velocity of sound is another vital input which goes into the hydrodynamics 
and control the evolution profile. Since electromagnetic probes are emitted throughout the spatial and temporal 
extent of the fireball, the spectra is in fact sensitive to each of these factors. We will investigate the dependence 
of the dilepton spectra in the low mass region on the equation of state (EoS) obtained from two different 
schemes. 

It is well known that the average magnitude of radial flow at the freeze-out surface can be extracted from the 
Pt spectra of various hadrons. However, hadrons being strongly interacting objects can bring the information 
of the state of the system when it is too dilute to support collectivity i.e. the parameters of collectivity 
extracted from the hadronic spectra are limited to the evolution stage where the collectivity ceases to exist. 
These collective parameters have hardly any information about the interior of the matter. On the other hand 
dileptons (and real photons) are produced and emitted from each space time point. Again, dileptons have an 
additional advantage because in this case we have two kinematic variables - of these, the pt spectra is affected 
by the flow in the system but the pt integrated invariant mass (M) spectra remains unaltered. Therefore, 
if one can identify a domain of M where QGP contributions dominate then the flow from the interior of the 
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QGP phase may be extracted by studying the My spectra at this M window i.e. a judicious choice of px and 
M windows will be very useful to characterise the flow in QGP and hadronic phases, which will shed light on 
the time evolution of the collectivity in the system [19]. The effective slope of the dilepton spectra have been 
extracted both for the RHIC and LHC energies to demonstrate these aspects. 

In this work we have thus focused on two kinds of many body effects which are inherent in the evaluation of 
dilepton spectra from heavy ion collisions. The one concerning the in-medium vector meson spectral function 
is microscopic in nature originating from the dynamics of effective hadronic interactions taken from chiral 
perturbation theory and the other is a macroscopic manifestation which results in collective flow of the fireball 
described by relativistic hydrodynamics. Both of these effects can be observed by studying the low mass lepton 
pairs with appropriate choice of invariant mass and transverse momentum windows. We demonstrate this by 
evaluating the lepton pair spectra for RHIC and LHC conditions. 

The organisation of the paper is as follows. In section 2 we shall discuss the details of calculation of the 
spectral function of vector mesons and subsequently the differential dilepton rate. Space-time evolution of the 
evolving matter will be discussed in section 3. Space-time integrated dilepton yields are presented in section 4 
and finally we conclude in section 5. 



2 Dilepton emission rate 

The rate of dilepton production rate from a thermal medium produced in heavy ion collisions is well known to 
be given by [20J 

° 2 T L (M 2 ) f BE (qo)g^ W, u (q ,q) (1) 



d*qd*x ~ 6rrV V "u» 

where the factor L(M 2 ) = (1 + 2mf/M 2 ) (1 — 4m 2 /M 2 ) 1 / 2 is of the order of unity for electrons, M(= ^fq 2 ) 
being the invariant mass of the pair and the electromagnetic (e.m.) current correlator W^ u is defined by 

W, u (Qo,q) = J d 4 xe^([j; m (x),Jr(0)]) (2) 

Here J^ m (x) is the electromagnetic current and () indicates ensemble average. For a deconfined thermal medium 
such as the QGP, Eq. ([1]) leads to the standard rate for lepton pair production from qq annihilation [21] at 
lowest order. 

At low values of the invariant mass M, the electromagnetic current of quarks J^ m = §W7 M « — \d"f^d — gS7/* s 
may be decomposed as 

JtT = J£ + J£/3-Jj/3 (3) 

where the vector currents 

J£ = ^(u-y^u + d-y^d) 

4 = W (4) 

are named by the lowest mass hadrons p°, u and <fi in the corresponding channel. Defining the correlator of 
these currents Wf^f^ analogously as in ([2]) we write, 

*V = W?, + W»J9 + W}J9 . (5) 
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Figure 1: The Dyson equation for the p propagator, the small circles indicating terms of higher order in the 
series. 

One now has to specify the coupling of the currents to the corresponding vector fields. For this purpose we 
write, in the narrow width approximation [22], 

(0\J e ™(0)\R) = F R m R e p (6) 

where R denotes the resonance in a particular channel and e M is the corresponding polarisation vector. The 
coupling constants Fr are obtained from the partial decay widths into e + e~ through the relation 

Fr = — w — (7) 

yielding Fr=0.156 GeV, 0.046 GeV and 0.079 GeV for p, uj and (ft respectively. Eq. (J6|) suggests the operator 
relations 

JP(x) = F p m p VP (x), JZ(x) = 3F u m u V£(x) etc. (8) 

where V/?^ (x) denotes the field operator for the p(io) meson. The vector correlator in the p channel for example, 
is then given by 

W>, = F'mlJd'xe^dV^VSiO)]) 

= K p p^(q ,q) (9) 

where pp\ u is the spectral function of the p meson in the thermal medium and K p = F p m 2 p . Defining the spectral 
functions of uj and <j> as in (0) , the dilepton rate is given by 

dN a 2 m 

^ = -^fBE(qo)g fMU [K p pP lv (qo, q) + K^ u {q , q) + K^qo, q)] (10) 

2.1 Spectral function in terms of self energy 

The p meson plays the most significant role in the low mass region. We thus begin with the complete p 
propagator which in the real time formulation of thermal field theory is a 2 x 2 matrix defined by 

D$(qo,q) = i J rf 4 xe^(r c y;(x)y I f(o)) a6 (11) 

where the superscripts (a, b = 1, 2) are thermal indices and T c denotes time ordering with respect to a contour 
in the complex time plane. In the case of a symmetric contour this matrix function factorises as [23J 

D >*- u [ o -d;J u > u -{ ^ vrr^i) (12) 
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and so is given essentially by a single analytic function, the diagonal component D^ u . This function admits a 
Kallen Lehmann representation given by |24J 

TSrM-rg , Mo ' ?> l , us) 

in terms of the spectral function defined in Eq. ([9"]). From this we obtain 

p£U9o, 9) = 2e(g )ImL>^(g , q) (14) 

Thus the dilepton rate is essentially given by the imaginary part of the p propagator which we will now elaborate 
in the following. 

The complete p propagator in the medium is obtained from the Dyson equation 

DfM = D$ ah {q) - Df x ac ( q W d ^(q)DfM (15) 

pictorially seen in Fig. [TJ D^} ab (q) is the free thermal propagator and IK* (g) denotes self energy insertions 
both of which are 2x2 matrices in this approach. Here we have considered one-loop diagrams consisting of a 
pion and another hadron h where h stands for ir, oj, a\ ,hi mesons restricting thus to the non-strange meson 
sector. These self energy graphs provide corrections to the p meson propagation in the medium and will be 
evaluated in the next section. 

One can get rid of the thermal indices on diagonalising the thermal propagator and self energy matrices so 
that Eq. (|15p now takes the form 

D,M = D^Aq)-D%n Xa ( q )D av { q ) (16) 

the bars denoting the diagonal components as discussed above. (q) turns out to be the vacuum propagator 
given by 

D%) = (-9,u + 2 ~* (17) 

The complete propagator D^ u {q) can be obtained in terms of transverse and longitudinal components by 
writing 

~D/iv = P^D t + QpvDi (18) 

and 

IV = PfJHt + Q^Jli (19) 

where the projection operators are given by 

Pfw = -9^u + -^rp - ^u^, Q^ v = -^j-u^Uj,, q 2 = (u ■ q) 2 - q 2 , (20) 

being four- velocity of the heat bath and u lL = — u ■ qq^jq 2 ■ Using these in Eq. (fT6|) we eventually arrive 
at the solutions 

D t (q) = ~* Di(q) = \- ~* (21) 

q z — m* — Iit{q) 9 9 - mf, - q z lli{q) 

neglecting the non-transverse piece in D^(q). These are used in Eqs. (fT4"|) and (fTU|) to arrive at the dilepton 
emission rate 

dN a 2 

fBE(qo)[K p A p (q ,q)+K iU A OJ (qo,q) + K lp A 4> (q ,q)}. (22) 



d A qd A x n 3 q 2 
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where e.g. A p (= -g^ Ppu (go, q)/G) is given by 



A. - - 

p 3 



+ 



2 £ Imllf 



{g 2 -mj-Yl Renf ) 2 + IndTf) 2 

g 2 imnf 



(g 2 -m 2 -g 2 Y^ Renf ) 2 + g 4 (£ Imnf) 



(23) 



the sum running over the ir — h loops. Thus, the dilepton emission rate in the present scenario boils down to the 
evaluation of the self energy graphs shown in Fig. [TJ The real and imaginary parts of the self energy function 
can be obtained from the 11-component as |23j 

Ren £L(<?0, g) = ReU^(g , g) = RelT^o, q) 

Imll^teo, g) = e(q )lmU(qo, g) = tanh( J Sg /2)Imn^(g , g) (24) 
where 11^ denotes the retarded self-energy. 

2.2 Evaluation of p self energy 
A. Hot meson gas 

The 11-component of the thermal self-energy matrix for the ir — h loops is given by 

n}i,(ff) = i J -^N^q,k)D^{k)D^{g-k) , (25) 
where D n (g) are the thermal propagators [23] and the factor N pu includes tensor structures associated with 



the two vertices and those in the vector propagator. In Ref. |14| . this was evaluated in detail using interaction 
vertices from chiral perturbation theory to obtain the imaginary part of the retarded p self energy. For positive 
values of go this is given by 



lmU p7T (g ,g) = -ir J 



d 3 k 

x 



(27r) 3 4w 7r w /i 

[N(k = £%)(! + n(oj n ) + n{io h ))5{go - u n - u h ) + N(k = -w 7r )(n(w 7r ) - n(uh))5(q + uj w - uj h )} .(26) 



where n is the Bose distribution function and the energy variables are co^ = y m 2 + k 2 , lo^ = y m 2 + (q — k) 2 . 
Here we have suppressed indices corresponding to the longitudinal and transverse components which have been 
evaluated as in [TJ]. We recall that the regions in the qo plane in which the two terms are non- vanishing give 
rise to cuts in the self energy function. These are controlled by the corresponding delta functions. The first 
term is non-vanishing for g 2 > (m^ + m^) 2 producing the unitary cut and the second term is non-vanishing for 
g 2 > (rrih — m^) 2 giving the Landau cut. The unitary cut is also present in vacuum but the Landau cut appears 
only in the medium and arises from scattering of p with the particles present there. With the aid of the delta 
functions the integration over k is easily performed. 

Note that the heavy mesons in |14j have been treated in the narrow width approximation. In the following 
we consider the conventional prescription (see e.g. |25j ) to take into account the decay widths of the h\ and a\ 
mesons without disturbing the analytic structure discussed above. Here, the self energy functions are convoluted 
with the (vacuum) spectral functions of the heavy mesons (h) as 
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(m h +2r h ) 2 




(27) 



h J 



{m h -2V h ) 2 



with 



(m h +2F h ) 2 




tt ml — m? + irhhT{rhh) 



(28) 



(m h -2T h ) 2 



where are the measured [26] decay widths and T(- v /i = rhh) are the corresponding ones calculated using 
hadronic interactions. As a consequence of this convolution the clear distinction between the regions of non-zero 
imaginary part will be smeared and for sufficiently large width might even appear continuous. However, for a 
particular value of the mass rhh, the imaginary part receives contribution from only one of the cuts. 

B. Effect of baryons 

It is known that the scattering of the p with baryons also contribute to the broadening of its spectral function 
even in the case of zero baryon density, and thereby increase the production of lepton pairs at invariant masses 
below the p peak. The relevant quantity in this case is the p self-energy due to baryon loops which has 
been evaluated at finite temperature and density using well established hadronic interactions in a many body 
approach [27] . 




The p self-energy at finite temperature and baryon density has also been estimated in terms of empirical 
scattering amplitudes |28j for p scattering from nucleons using experimental data. The self energy of the p is 
obtained from these amplitudes using a low density approximation as 



The scattering amplitude fp'™'(s) in the low energy region is described in terms of N* and A baryon resonances 
and a Regge model is employed at higher energies |28j . Since the spectral function in this approach seems to 
agree with the many body approach as shown in [29] we will adopt this formalism to estimate the baryonic 
contribution to the self energy. 

The total p self energy in hadronic matter containing mesons and baryons is then written as 



The spectral function as a function of M is plotted in Fig. [2] for a range of temperatures at a baryonic 
chemical potential ps = 30 MeV relevant for RHIC energies. We observe significant broadening at higher 
temperatures with almost negligible shift in the pole mass. 

2.3 The dilepton emission rate at temperature T 

We next plot in Fig. [3l upper panel, the dilepton emission rate keeping only the p contribution in Eq. (|22p in 
which we show, for illustrative purposes, the relative contributions from the cuts in the tt — h loops keeping 
only one of them at a time. The unitary and Landau cuts for the tt, u, hi and a± are seen to contribute with 
different magnitudes for different values of the energy and three momenta of the off-shell p. In the time-like 
region, in the vicinity of the (bare) rho mass the imaginary part of the self energy from a particular loop 
receives dominant contribution from only one of the cuts. The tt — tt loop for example, has only the unitary 
cut and this contributes most significantly to dilepton emission near the p pole. In contrast, the Landau cut 
contribution from the tt — u loop is dominant up to about 400 MeV. Since this cut ends at M = — m n and 
the unitary cut starts at M = + m n there is no contribution at the p pole. However, the unitary cut of the 





(29) 



(30) 
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Figure 2: The p spectral function averaged over the transverse and longitudinal polarisations for <f=300 MeV 
and p B = 30 MeV. 
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M (GeV) 

Figure 3: Upper panel shows contributions from the discontinuities of the self-energy graphs to the dilepton 
emission rate at T = 175 MeV and ps = 30 MeV. L and U denote the Landau and unitary cut contribution. 
Lower panel shows contributions from the mesons and baryons. 



lott loop could make a significant contribution as seen in |141 130] . The Landau cut for the n — a\ self-energy 
extends up to about 1100 MeV and makes a substantial contribution both at and below the p pole. The unitary 
cut starts at a much higher value of M and hence does not make a significant contribution to the p spectral 
function. We also show the effect of convolution over the width of the a\ as discussed above. As expected, 
the contributions from the Landau and unitary cuts are now joined by a continuous line, the boundaries being 
smeared out due to the substantial width of the a\. While analysing the different contributions one must keep 
in mind that the total contribution from the different loops to the spectral function is not a linear sum of the 
individual contributions as is clear from the definition given in Eq. (|23p . This is seen in the lower panel where 
the cumulative contribution to the lepton pair yield is shown for the it — ir and ir — h loops. Also shown is the 
enhancement in yield obtained by including baryons at ps = 30 MeV for RHIC energies. 

We have also included dilepton emission from the uj and the <f>. The width of the to at finite temperature is 
taken from the calculation of Ref [31] where a framework similar to the one employed here has been used. For 
the 4> only the vacuum width has been considered. In addition, the spectral function of the p and uj has been 
augmented with a continuum contribution given by 

3qK - E *T PU 

V=p,uj 

where the continuum part is parametrised as [16\ [22] 

1 (. ats\ 1 



K = i- 1 + — T~. 7 ttt (32) 



8tt V 7r / 1 + exp(wo — qo)/S 



with coq = 1.3 (1.1) for p (w) GeV and 5 = 0.2 for for both p and uj. The continuum contribution for the uj 
contains an additional factor of 1/9. The dilepton emission rates from QGP and hadronic matter have been 
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plotted in Fig.[4]at a temperature of 175 MeV and baryonic chemical potential of 30 MeV. We observe significant 
enhancement in the dilepton yield in the mass region below the p pole compared to vacuum. This rate has been 
used in the analysis of the dimuon spectra obtained from In-In collisions at 17.3 GeV at CERN SPS [32J. The 
calculations show a reasonable agreement with the invariant mass spectra for different pr ranges as well as the 
Mt spectra for different M bins. 

3 Space time evolution 

Thus far we have discussed the dilepton emission rates for a given temperature. In a HIC the temperature 
corresponds to that of local equilibrium and is hence a function of position and time. For a quantitative 
evaluation of the dilepton yield one has to convolute the static rate over space and time resulting in an integrated 
yield which is a superposition of contributions from a range of temperatures. Ideal relativistic hydrodynamics 
is used to carry out this scheme in which a number of inputs are required. These are discussed in the following. 

3.1 Equation of state and initial conditions 

The evolution of the fluid is governed by the energy momentum conservation equation 

<9 M T^ u = (33) 

where T^ v = (e + P)u^u u + g^ u P is the energy momentum tensor for ideal fluid. This together with a relation 
connecting the pressure P and the energy density e, known as the EoS provides a closed set of equations [33J. 
With cylindrical symmetry and longitudinal boost invariance [34] . these set of equations provide the energy 
density and transverse velocity as a function of the proper time and the radial coordinate. The initial energy 
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density and radial velocity profiles which go as inputs are taken as: 



C fo' r > = l + e (r-flA)/f (34) 

and 

v(n,r) = (35) 

5 (~ 0.5 fm here) is a parameter, known as the surface thickness. The static rate, dR/dM 2 at M = 0.4 GeV 
with hb = 30 MeV and (jlb = differs by less than 10%. The difference will be even smaller at LHC energies 
because of the smaller value of p,B at central rapidity region. Therefore, the effects of hb on space-time evolution 
is neglected in the present work. 

It is natural to expect that different EoS's will govern the hydrodynamic flow quite differently [T7] and 
as far as the search for QGP is concerned, the goal is to look for distinctions in the observables due to the 
different EoS's (corresponding to the novel state of QGP vis-a-vis that for the usual hadronic matter). It is 
thus imperative to understand in what respects the two EoS's differ and how they affect the evolution in space 
and time. In order to check sensitivity with the equation of state (EoS), we have considered two scenarios: (a) 
hadronic resonance gas (HRG) with all hadrons up to mass 2.5 GeV for the hadronic phase along with a bag 
model EoS for the QGP phase and (b) EoS obtained from lattice QCD calculations (LQCD) [35J. 

One of the most important parameters that go into the space time evolution are the values of the initial 
temperature and the thermalisation time. There are indications that QGP thermalises quite early at RHIC 
energies. In case of isentropic expansion the experimentally measured hadron multiplicity can be related to the 
initial temperature and thermalisation time by the following equation |36| : 



\ 2-7T 4 ,dN .. . . 

T '^ )T ' - m¥w^ { ^ (bm)) (36) 

where {dN/dy(b m )) is the hadron (predominantly pion) multiplicity for a given centrality class with maximum 
impact parameter b m . R± is the transverse dimension of the system, T{ is the initial thermalisation time, ("(3) 
is the Riemann zeta function and = (ir 2 /90) is related to the degeneracy (<?&) of the system created. The 
hadron multiplicity resulting from Au + Au collisions is related to that from pp collisions at a given impact 
parameter and collision energy by 

<^(U> = [(1 - x)(N part (b m ))/2 + x(N coll (b m ))} ^ (37) 

where x is the fraction of hard collisions, (N par t) and {N co u) are the average numbers of participants and collisions 
respectively evaluated by using Glauber model. dN^/dy = 2.5 - 0.25Zn(s) + 0.023/n 2 s, is the multiplicity of 
the produced hadrons in pp collisions at centre of mass energy, -y/i [37]. Assuming 10% hard (i.e. x = 0.10) and 
90% soft collisions for initial entropy production the value of dN^/dy turns out to be about 2.43 at = 200 
GeV. For RHIC energy, we take T{ = 320 MeV with initial time Tj = 0.2 fm/c which acts as inputs to the 
hydrodynamic evolution. 

At LHC the measured values of dN^/dy for y/s^ = 900 GeV, 2.36 TeV and 7 TeV are 3.02, 3.77 and 6.01 
respectively [38]. The value dN^/dy at y/s^N = 5.25 TeV is obtained by interpolating the above experimental 
data mentioned above. Assuming x = 0.2 in Eq. (I37j) we obtain dN/dy = 2607 in Pb+Pb collision for 0-10% 
centrality. For r, = 0.1 fm/c we get Tj = 686 MeV. 

For the space-time picture, we thus work in the following scenario. An equilibrated QGP is formed at 
initial temperature (time) Tj(rj), the system then cools due to expansion and when the temperature reaches 
T c it undergoes a phase transition from QGP to hadrons. After the completion of the phase transition the 
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hadronic matter cools and eventually freezes out first chemically at a temperature T c ^ and then kinetically at a 
temperature Tp. The transition temperature is taken as T c ~ 175 MeV. The other inputs which goes into the 
calculations are chemical (T c ^) and kinetic freeze-out (Tp) temperatures. The kinetic freeze-out in the system 
occurs when both the elastic and in-elastic collisions stop i.e. the freeze-out takes place when the collectivity 
in the system ceases to exist. The value of Tp can be constrained from the hadronic pp spectra pQ. In the 
present work we take Tp = 120 MeV which reproduces the pp spectra of pions, kaons reasonably well [39J. The 
ratios of various hadrons measured experimentally at different ^/snn indicate that the system formed in heavy 
ion collisions chemically decouple at T c ^ which is higher than Tp [3D]. Therefore, the system remains out of 
chemical equilibrium from T c h to Tp. The deviation of the system from the chemical equilibrium is taken in 
to account by introducing chemical potential for each hadronic species [41 . The chemical non-equilibration 
affects the yields through the phase space factors of the hadrons which in turn affects the productions of the 
EM probes. The chemical potential, fij for the hadronic species j as a function of T have been taken from 
Ref. 02]: 

n j( T ^j) rijjTch,^ = 0) 

a(T,{ Mi }) s(T ch ,{ f x j } = 0) {M) 

where rij is the density of hadron j contains direct as well as contributions from resonance decays. The fij 
is a function of T and it vanishes at T = T c h (= 170 MeV here). Therefore, the space time evolution of \ij 
is dictated by the evolution of T. The chemical potentials of pions, uj, hi, a\, (j) and proton enters through 
their thermal distributions as a fugacity factor. The values of the respective chemical potential at the kinetic 
freeze-out temperature, Tp = 120 MeV are fi n = 68 MeV, ^ = 179 MeV, fi hl = 204 MeV n ai = 204 MeV, 
= 252 MeV /i proton = 258 MeV. 

We have also included the contribution to the dilepton yield from the decays of p mesons at freeze-out using 
the Cooper-Frye formula [43j as follows. For a special case of unstable vector mesons we need to know the 
thermal phase space factor corresponding to an unstable Boson which is given by 



where E = \/p 2 + M 2 , g is the statistical degeneracy and A p (M) is the (vacuum) spectral function of the vector 
meson under consideration. For stable particle A p {M) reduces to a Dirac delta function and consequently the 
usual phase space factor for a stable particle is recovered upon integration over dM 2 . Therefore, the rap 
distribution of dimuons from vector meson decay after the freeze-out is given by [8] (see also [44]): 



M T dM T dM 2 dy 



2ir j dr j dr/ j d(f>rr 

( dr 
x I Mrcosh(y — rf) — —ppcoscj) 

x j unstable^ V->-/i+^- /Ty*' 

(40) 

where Ty* is the total decay width of the vector meson V. 

4 Dilepton spectra at RHIC and LHC 

We begin by plotting the space-time integrated invariant mass spectra of dileptons. In Fig. Owe plot the yield 
of lepton pairs from the hadronic matter (HM) containing the effects of chemical non-equilibration both in 
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Figure 5: Invariant mass distribution of dileptons from hadronic matter (HM) for modified and unmodified p 
meson. Scenario (a) has been used for EoS. 

the pole as well as in the continuum part of the spectral function, evaluated with and without the modified 
p spectral function discussed in Section 2 for RHIC. The enhancement in the region 0.1 < M < 0.7 GeV is 
purely a medium effect and is a contribution from the Landau cut of the tt — lj, h\, a\ loops. In simple terms, 
the yield in this region results from the scattering of the p with thermal gj, hi, a\ and pions |14j . As seen 
from Eq. (|26p these contributions are weighted by the Bose distribution functions for the pions (minus a much 
smaller contribution from the heavy mesons). In contrast, the vacuum spectral function naturally starts from 
the 2m n threshold coming from the unity in the unitary cut contribution. The (small) kink at 0.42 GeV in 
this curve is due to the 3m n threshold for uj production. The enhancement in the yield due to medium effects 
is ~ 20 for M around 400 MeV. Next we show in Fig. [6] the dependence of the yield from the two phases 
on the EoS. Dilepton radiation from hadronic phase outshines the emission from quark matter for M up to 
<f> mass. Since we have included the continuum in the p and u channels we ignore the four pion annihilation 
process [45] to avoid double counting. The contributions from quark matter phase dominates over its hadronic 
counter part for both the EoS for M beyond </>-peak. This fact may be used to extract various properties i.e. 
average flow, temperature etc. of quark matter and hadronic matter by selecting M windows judiciously. The 
dilepton yield from hadronic matter is observed to be larger when the HRG EoS is employed in comparison to 
LQCD. This can be understood in terms of the velocity of sound Vg(= dP/ale evaluated at constant entropy) 
which controls the rate of expansion. For EoS of the type (a) Vg ~ 1/3 in the QGP phase which is larger than 
the value of the corresponding quantity for EoS of the type (b). Therefore, the rate of expansion in the scenario 
(b) is comparatively slower, allowing the QGP to emit lepton pairs for a longer time resulting in greater yield 
for LQCD EoS. In contrast, for the EoS (a), the lower value of v s for the hadronic phase results in a slower 
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Figure 6: The invariant mass distribution of lepton pairs for RHIC initial condition with HRG and LQCD EoS. 



cooling and hence a larger yield. Also shown for comparison is the yield from the decays of p mesons at the 
freeze-out for the two types of EoS used. The yield from this source is much smaller and we will not consider 
it any further. 

Since the M spectra is invariant under flow we now turn to the My(= \Jpt + ^aL) spectra to study this 
aspect. Fig. [7] shows the My spectra of lepton pairs at RHIC energies. Here the differential yield is integrated 
over small bins of the pair invariant mass (from M\ to M2) and plotted against My — M av which is actually a 
measure of the average kinetic energy (KE) of the pair, M av being the average mass (= [M\ + M<2\/2) of the 
bin. The average value of Mt for a static system at a temperature T is given by (Mt) ~ M + T. Therefore, 
the average KE ~ T, is the slope of the Mt distribution. Initially, the entire energy of the system formed 
in HIC is thermal in nature and with progress of time some part of the thermal energy gets converted to the 
collective (flow) energy. In other words, during the expansion stage the total energy of the system is shared by 
the thermal as well as the collective degrees of freedom. As a consequence, unlike the invariant mass spectra, 
the Mt (or pr) spectra is heavily influenced by the collective flow and the average KE or the inverse slope may 
be written as T e ff = T + l/2M av v T , where vt is the average radial flow velocity. The Mt spectra of dileptons 
for various M-bins have an exponential nature, the inverse slope providing an effective temperature, T e ff. It 
is important to mention at this point that for a radially expanding system the T e ff has an explicit (linear) M 
dependence as mentioned above. However, it has also an implicit M dependence even when vt = (i.e. with 
longitudinal expansion only) because it is expected that the high (low) M pairs predominantly emit from the 
high (low) temperature or early (late) time zone. For a radially expanding system the M dependence of the 
T e ff is stronger than for a system which expands longitudinally only. 

In Fig. [8] we have plotted the effective temperature versus M av for various mass windows of the lepton 
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pairs at RHIC energies, evaluated with the in- medium spectral function of the vector mesons. Also shown by 
a filled square is the value of T e ff for the vacuum case in the window 0.4 < M < 0.6 where there is substantial 
difference between the yields in free and medium cases as seen in Fig [5j The slope of these curves measure the 
average temperature and the flow of the matter. 

Let us try to understand the non-monotonic variation of the inverse slope with M av depicted in Fig. El 
In Fig. [6] it is shown that the high M (above <fi peak) pairs originate predominantly from the partonic source 
and the low M (below p mass) domain, although outshine by the radiation from hadronic source, contains 
non-negligible contributions from quark matter i.e. the low M region contains contributions both from the 
hadronic as well as QGP phases. Now, the collectivity (or flow) in the system does not develop fully in the 
QGP because of the small life time of this phase which means that the radial velocity extracted from the high 
M region is small. Here the temperature decreases mainly due to longitudinal expansion and consequently, the 
effective slope decreases slowly with decreasing M av . In contrast, the lepton pairs with mass around p-peak 
dominantly originate from the hadronic source (which appears in the late stage of the evolving system) and are 
significantly affected by the flow resulting in higher values of vt and hence a higher T e ff. At still lower values 
of M(< m p ) the contributions from hadronic phase is considerable. Therefore, the Mt spectra in this domain 
is affected by the flow significantly i.e. the value of vt is large but can not be as large as in the p-peak region 
despite the substantial medium induced enhancement of the hadronic sources since this domain also contains 
contribution from the QGP. But for low M(< m„) the T e ff is smaller compared to the p peak region because 
of the smaller values M. The value of T e ff increases linearly with M up to the p peak, a behaviour typical 
of inverse slope extracted from the transverse momentum spectra of hadrons with different masses. Thus the 
values of T e g for M below and above the p-peak are smaller compared to the values around the p peak even 
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Figure 9: Dilepton yields for HRG EoS and LQCD EoS. The initial condition is taken for LHC energy. 

in the presence of medium effects, resulting in the non-monotonic behaviour as displayed in Fig. [8j for 0.5 < 
M(GeV) < 1.3. 

The slope of the Mt spectra is connected with the average collective flow. It is well known that the average 
magnitude of radial flow at the freeze-out surface can be extracted from the pr spectra of the hadrons. However, 
hadrons being strongly interacting objects can bring the information of the state of the system when it is too 
dilute to support collectivity i.e. the parameters of collectivity extracted from the hadronic spectra are limited 
to the evolution stage where the collectivity ceases to exist. These collective parameters have hardly any 
information about the interior of the matter. On the other hand the dileptons are produced and emitted from 
all space time points. Therefore, the value of vt estimated from the dilepton spectra will be lower than the 
value extracted from the hadronic spectra Indeed, the values of vt estimated from the slopes of the curve 
is 0.25 for the M domains 0.5 < M (GeV)< 0.77. This value is much smaller than the value of vt extracted 
from the hadronic spectra [IT]- The dashed line in Fig. [8jis obtained by setting vt = 0. The results indicate 
that the observed (solid line) rise (for 0.5 < M(GeV)< 0.77) and fall (for 0.77 < M (GeV)< 1.3) are due to 
radial expansion of the system. However, the rise in large M domain is due to cooling of the system due to 
longitudinal expansion - which is described as the implicit M dependence of T e ff above. 

The invariant mass spectra of lepton pairs is displayed for LHC initial conditions in Fig. [9j Although, the 
results are qualitatively similar to RHIC, quantitatively the yield at LHC is larger. This is because of the larger 
four- volume of the system to be realised at LHC resulting from a higher value of Tj for fixed T c and Tp. Similar 
kind of enhancement is also reflected in the transverse mass distributions of the lepton pairs at LHC (Fig. ITU]) . 

Finally the variation of inverse slope of the Mt distributions with M av for LHC is depicted in Fig. [TTJ 
The values of T e ff for various M-bins are larger than RHIC because of the combined effects of large initial 
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temperature and flow. In fact the value of vt for 0.5 < M(GeV)< 0.77 is ~ 0.52 compared to 0.25 at RHIC. 
The radial flow in the system is responsible for the rise and fall of T e ff with M av (solid line) in the mass 
region (0.5 < M(GeV)< 1.3), for vt = (dashed line) a completely different behaviour is obtained. This type 
of non-monotonic variation of T e // can n °t be obtained with a single dilepton source [48]. Therefore, such 
non-monotonic variation of the inverse slope deduced from the transverse mass distribution of lepton pairs with 
average invariant mass is an indication of the presence of two different phases during the evolution of the system. 
Thus, such variation may be treated as a signal of QGP formation in heavy ion collisions. A comment on the 
thermal emission rate from the QGP is in order here. We have considered the lowest order processes to evaluate 
the the lepton pair productions form QGP. However, we have checked that the inclusions of the processes of 
order 0(aa s ) \49 \ \50 \ ISTj changes the slope parameter, T e ff at low M by negligible amount. It is less than 1% 
for M ~ 0.3 GeV and at higher M it is vanishingly small. 

The other sources of dileptons e.g. from the Drell-Yan (DY) mechanism and charm decays may provide 
significant "background" to the thermal productions at high mass region (2 < M (GeV) < 6 [52]) which are 
neglected here because in the present work we focus mainly on the low mass regions. Moreover, the contributions 
from the DY process and charm decays from proton+proton (pp) collisions may be used to estimate the similar 
contributions from heavy ion collisions at the same colliding energy by appropriately scaling pp data by the 
effective number of nucleon+nucleon collisions in nuclear interaction. The contributions from the p at the 
freeze-out surface has been evaluated and it is found to be small. 

5 Conclusion 

In this work we have attempted to bring out distinguishing features stemming from many body effects in the 
lepton pair yield from relativistic heavy ion collisions. On the microscopic side we have used a p spectral function 
evaluated at finite temperature in the framework of real time formalism of thermal field theory using interaction 
vertices from chiral perturbation theory. The effect of baryons have also been included in the empirical approach 
of Ref. [28] using resonance dominance in the forward scattering amplitude. A significant enhancement of the 
dilepton yield in the region below the nominal p peak has been obtained from the combined effect of mesons and 
baryons. The space-time evolution of the dilepton emission rate using relativistic hydrodynamics with chemical 
freeze-out reflects this enhancement in the invariant mass spectra. 

On the macroscopic side, it is argued that the non-monotonic variation of the inverse slope deduced from the 
transverse mass distribution of lepton pairs for various values of the average invariant mass is an indication of 
the presence of two different phases during the evolution of the system. Thus, such a variation may be treated 
as a signal of QGP formation in heavy ion collisions. 
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